Effective anisotropy of thin nanomagnets: beyond the surface anisotropy approach 
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We study the effective anisotropy induced in thin nanomagnets by the nonlocal demagnetization 
field (dipole-dipole interaction). Assuming a magnetization independent of the thickness coordinate, 
we reduce the energy to an inhomogeneneous onsite anisotropy. Vortex solutions exist and are 
ground states for this model. We illustrate our approach for a disk and a square geometry. In 
particular, we obtain good agreement between spin-lattice simulations with this effective anisotropy 
and micromagnetic simulations. 
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PACS numbers; 75.10.Hk, 75.70.Ak, 75.40.Mg, 05.45.-a 



I. INTRODUCTION 

Magnetic nanoparticlcs and structures have recently 
attracted a growing interest for their phvaical properties 
and a number of possible applications E'B'El For example 
the vortex (ground) state of a disk-shaped nanoparticle 
could provi|de high density storage and high speed mag- 
netic RAMq. The theoretical models for these systems 
have been known for some timeBu and include the nonlo- 
cal demagnetization field. At microscopic level this field 
is due to the dipolar interaction 



n.m 



(1) 



Here Sn = {S!^,S^,S^) is a classical spin vector with 



fixed length S on the site n 



,nz) of a three- 



dimensional lattice. The summation runs over all mag- 
nets {n,m), and Tmn = fn — I'm- The parameter 
D = n\g^ is the strength of the long range dipolar inter- 
action and g is the Lande-factor. 

In the past analytical studies have been mainly limited 
to assuming a homogeneous demagnetization field dis- 
tribution, uniform Stoner-Wohlfarth thcoryQ and near- 
uniform Brown's linear analysisEI. Recent advances in 
nanotechnology and computing power established the 
complexity of magnetization distribution in nanoparti- 
cles. For example square nanoparticles exhibit buckliaft 
states, flower states, apple states, leaf states etcJ3tl't2l 
when their size exceeds the siiiftle.-domain limit. In disk- 
shaped pacticies vortex stateaJ'EJ, edge fractional vortex 
states etdlaGJ appear. Some of these complex states 
can be obtained by a small perturbation of a homp^ 
geneous state. For example Cowburn and WellandEdl 
showed that dipolar interactions cause flower and leaf 
states in square nauioparticles, which was confirmed by 
direct experiments^. However the linear analysis does 
not work for topologically nontrivial states like kinks, 
vortices etc. One possibility to study these structures in 



nanomagnets is the Ritz variational method. It was ap- 
plied to japalyze the vortex structure of the disk-shaped 
nanodotElEil. A disadvantage of this method is to limit 
the solution to a certain class of minimizers, so that one 
can usually study only one type of excitation. Linear 
waves are left out together with their coupling to the 
main excitation. 

The various regimes were studied in Refs. [l5| , [l6| , p7|jl^ , 
p^ , pO 21, The import ant length scale is the magnetic ex- 
change length t = \/ k/ AiiMg where A is the exchange 
constant and Mg is the saturation magnetization. De- 
pending on the relation between the film diameter 2i?, its 
thickness h and € many scaling limits can be analyzed, 
see Ref. |l^ for an overview. Probablji. the first rigorous 
study was made by Gioia and Jamesll3 who showed that 
for an infinitesimally thin-film {h/R ^ 0, — > const) 
the magnetostatic energy tends to an effective 2D easy- 
plane anisotropy energy. In this case the ground state 
is a homogeneous in-plane magnetization statcEj. This 
effective easy-plane anisotropy has a simple magneto- 
static interpretation. The sources of magnetostatic field 
are volume and surface magnetostatic charges. For thin 
structures one can neglect the volume charges. Face sur- 
face charges contribute to the energy density as 2ttM'^ 
which is the same term oiie would get with an effec- 
tive easy-plane anisotropycj. In the case h/R 1 and 
£^ ^ 2/i i? I ln(/i/2i?)| the magnctizaliticL. develops edge 
defects, including fractional vortices Ejcj'EJ This prob- 
lem has a boundary constraint and an interior penalty. 
It is relevant for typical Permalloy (Ni8oFe2o, Py) disks 
where we have h 20nm, 2R ^ lOOnm and £ 5.3nm. 

It was shown in Refs. |lj,|l^ that in the limit h/R 
under the scaling 
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the full 3D micromagnetic problem reduces to a much 
simpler 2D variational problem where the magnetostatic 



2 



energy tends to the effective surface anisotropy term 



£;surf = J{S-TfdS 



(3) 



on 



where r is the local tangent vector on the domain bound- 
ary dn. In this case the magnetization S has no out of 
plane component {Sz = 0) and does not develop walls 
and vortices. 

To study nanomagnets with curling ground states, here 
we develop a new analytical approach. We split the 
dipole-dipole spin interaction into two parts. The 
first one is an on-site anisotropy with spatially depen- 
dent anisotropy coefficients. The second part represents 
an effective dispersive interaction. The anisotropy inter- 
action consists of two terms: an easy-plane anisotropy 
and an in-plane anisotropy. We show that the vortex 
state minimizes the effective in-plane anisotropy. We also 
show that for ultra-thin nanomagnets {h/ 0) the in- 
plane anisotropy term reduces to the surface anisotropy 
(^. For the nonhomogeneous state which is our main 
interest, our approach is valid if 



i? > /i and i? > 



(4) 



In Sec. we introduce our discrete model together 
with the dipolar energy and adapt it to the plain-parallel 
spin-field distribution, which is our main simplification. 
We further simplify the model by considering only the 
local part of the dipolar energy, which results in an effec- 
tive anisotropy. In the continuum approximation of the 
system we get a local energy functional with nonhomo- 
geneous anisotropy coefficients, see Sec. II A. The disper- 



sive interaction is discussed in Sec. [IB. To illustrate our 



method of effective anisotropy we consider in Sec. Ill the 



disk-shape nanoparticle and study its ground state spin 
distribution. Our simple m odel d escribes exactly the ho- 
mogeneous state (see Sec. Ill A) and very precisely the 
vortex state (see Sec. IIIB ). In Sec. IV we confirm our 
analysis by numerical simulations. These are done first 



for the disk-shaped nanoparticle (Sec. IV A ) and then for 
the prism-shaped one (Sec. fV B[ ). We discuss our results 
in Sec. 



II. MODEL. EFFECTIVE ANISOTROPY 

We consider a ferromagnetic system described by the 
classical Heisenberg isotropic exchange Hamiltonian 



J 



(5) 



where the exchange integral J > and the summa- 
tion runs over nearest-neighbor pairs (n,n'). The total 
Hamiltonian is a sum of the exchange energy (|^) and the 
dipolar one (|l|). 

Our main assumption is that depends only on the 
X and y coordinates. Such a plane-parallel spin distribu- 
tion is adequate for thin films with a constant thickness 



h = NzttQ, (oo being the lattice constant) and nanopar- 
ticles with small aspect ratio. The exchange interaction 
can be written as the sum of an intra-plane 3^™'"''' term 
and an inter-plane one 'K]^°'^ 



riru 



1)J 



(6) 



Here and below the Greek index u = {rix, riy) corresponds 
to the XY components of the vector n ~ {nx,ny,nz). 
One can see that the inter-plane interaction is equivalent 
to an on-site anisotropy. The inter-exchange term gives 
a constant contribution, so it can be omitted. 

Let us consider the dipolar energy. Using the above 
mentioned assumption about the plane-parallel spin dis- 
tribution, the dipolar Hamiltonian can be written as (see 
Appendix ^ for the details): 

•^t^ (7a) 

Here the sum runs only over the 2D lattice XY. All the 
information about the original 3D structure of our system 
is in the coefficients A^^, Bfj^, and C^i/, 
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To gain insight into the anisotropic properties of the 
system we represent the dipolar energy (Wm as a sum 



rw-loc 



where 
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B^[{s^f-{syf]+2c^s^sy 



(9) 



is an effective on-site anisotropic energy and 



(10) 
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FIG. 1: (Color online) Arrangement of coordinates in 
the local reference frame. 



is the dispersive part of the dipolar interaction. Here we 
introduce the coefficients of effective anisotropy 



(11) 



The dipolar energy IK^°'^ contains only local interaction; 
it has a form of the anisotropy energy with nonhomo- 
geneous A^, Bi,, d,- In next sections we discuss these 
quantities. For this end we need to obtain the continuum 
limit of our model. 



A. Continuum description 

The continuum description of the system is based on 
smoothing the lattice model, using the normalized mag- 
netization 



m{r) = 



alMs 



r -Vr 



(12) 



where Ms is the saturation magnetization. The exchange 
energy, the continuum version of (^) is 



£cx = + ao) J d^x {Vmf , 



(13) 



where A = JMga^/D is the exchange constant. 

Now let us consider the dipolar energy and use its ap- 
proximate Hamiltonian (^). To present this energy in a 
standard phenomenological form one needs to transform 
the summation over the lattice to an integration over the 
volume. There is a singularity for 0. Using a 

regularization similar to the one in Ref. |2^, we find (see 
Appendix p| for details) that the local part of the dipolar 



energy is 



£<j nM'^h J d^xi^A{x,y) [l - Scos^ 9] 



(14) 



' si-D^e Re 



where we used the angular parameterization for the mag- 



netization: m 



cos 6 and m'^ 



sin^e"^. Here 



and below we dropped the loc superscript. One can see 
that the original nonlocal dipolar interaction results in 
an effective local anisotropy energy. The coefficients of 
effective anisotropy A and ® are nonhomogeneous: 
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S(x,y) 
J(P, h) 



APh Ah{P^ + h^)3/2 

l-l%{P,h)e-^'"da, 

P - + 



(15a) 
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(15b) 



where the Heavisidc function 0+ (x) takes the unit values 
for any positive x and zero values for x < 0. In Eqs. ( p^ ) 
we used the local reference frame 

x' = X + pcos{x + a), y' = y + psm{x + a), (16) 

which is centered at {x,y). The term P is the distance 
from this point to the border of the system, it depends 
on the azimuthal angle a and position {x, y), see Fig. ^. 

In the limit case of the pure 2D system (monolayer 
with /i = 0) the total energy, normalized by the 2D area 
§, takes a form 



W 



h=0 _ 



A/|§ao 
27r£2 



d ^ ^h=0 



w. 



h=0 



J d^x [{Vef + sin^ e{V(j)f] 
■t' = lld'x!^A'^='{x,y) [l-3cos^0] 

S'^=0(x,y)e2'(*->^)l|, 



(17) 



+ sin^ 9 Re 
A'^=°{x,y) 
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Here the exchange length ^ has the standard formB 



47rM2 



ao 



AttD 



(18) 



induces magnetic anisotropy was 
•2^' for a pure 2D spin system from a 



Note that the dipo" 
considered by Levy 
Taylor's series expansion of the spin field. 

The above case (^7|) has rather an academic inter- 
est. Below in the paper we consider another limit, when 
h ^ uq. In that case one can neglect the energy of the 
monolayer The total energy, normalized by the 

volume of the magnet, takes a form: 



W 



dfx 



(19a) 
(19b) 



Wd = ^ J d^xi^A{x,y) [1 - 3 0032 6*] 



sin2 6'Re 'B{x,y)e 



o2j(0-x) 



(19c) 



The effective anisotropy constants can be expressed as 
follows: 



A{x,y)^-J^ da 



(20a) 
(20b) 



nPM^ , -21n ^^ \ (20c) 

h F 

Let us discuss the magnetization distribution of the 
nanoparticle on a large scale. The equilibrium mag- 
netization configuration is mainly determined by the 
dipolar interaction, which takes the form of an effective 
anisotropy (19c). The coefficient A determines the uni- 
axial anisotropy along the z-axis. For thin nanoparticle 
this coefficient is always negative (with A —2/3 when 
h ^ Q) favoring an easy-plane magnetization distribu- 
tion in agreement with the rigorous calculationalj. The 
coefficient 25 is responsible for the in-plane anisotropy 
in the XY-plane. Assume that all spins lie in the plane 
corresponding to the thin limit case. The preferable mag- 
netization distribution in the XY-plane is the function 



minimizing the expression Re 



3g2j(</.- 



= X + 2 - 2 



(|l|). This 



(21) 



The angle (|2l] ) determines the in-plane effective 
anisotropy direction observed on a large scale, without 
exchange interaction and effective uniaxial anisotropy. 
The analysis of the S-term shows that the effective 



anisotropy favors such an in-plane spin distribution, al- 
ways directed tangentially to the border near the sample 
edge (see Appendix ^ for the details). This statemcajt 
agrees with results for the pure surface anisotropyllj. 
Finer details depend on the geometry of the particle so 
we need to distinguish the disk shape from the square 
shape. 



B. Dispersive part of the dipolar interaction 

In the continuum description (^) the dispersive part 
of the dipolar interaction (R^ takes the form 



- 3 K(r) - m'{r')f^ + B{r - r'){[m^(r) - m^(r')]^ 

- [m^(r) - my{r')f^ + 2C(r - r') [m^(r) - m^(r')] 
X [TO^(r) - m^(r')] . 

By applying the Fourier-transform 

rn(r-) = ^ / d''7m(q)e^''-, (22) 

and neglecting finite-size effects, the normalized dis- 
persive part of the dipole-dipole interaction AWd = 
A£d/(M|§/i) can be represented in the form 



\z\2 



\q ■ m„ 



(23) 



Here q = {qx,qy) is the two-dimensional wave vector, 
rn{q) is the Fourier-component of the two-dimensional 
magnetization m(r), and the function ©(q) is defined by 
the expression 



6(g) 



qh-l 



-qh 



qh 



(24) 



Note that Eq. (|2^) is obtained under assumption that the 
ortho- normalization relation 



1 



(27r)2 



d X e 



l{q~q')-r 



S{q - q') 



takes place. Being exact for the infinite domain, this re- 
lation is only approximate for the finite-size system. For 
qh ^ the function (|2|) takes the form (5(q) w qh/2. 
Therefore we expect our approach to yield the correct re- 
sults for the homogeneous and for weakly inhomogeneous 
states. For the general nonhomogeneous spin distribu- 
tion the effective anisotropy approach gives only approxi- 
mate results. In Sec. [II we verify our effective anisotropy 
model for disk-shapes nanoparticles. 
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FIG. 2: (Color online) Spacial dependence of the effective anisotropy constants yi [see Eq. (B7)] and 23 [sec Eq. (B9) 



III. DISK-SHAPE NANOPARTICLE 

Let us consider a cylindric nanoparticle of top surface 
radius R and thickness h. We introduce £ = h/{2R) the 
aspect ratio. Let us calculate first the effective anisotropy 
coefficients A and 53. For the circular system the coef- 
ficients A and ® depend only on the relative distance 
^. We calculated analytically the coefficients A and 23 
(see AppendixH) and these are presented in Fig. and 
Eqs. (p7|), (|B9|). First note that when e > 1 both 
anisotropy constants asymptotically do not depend on 

A{Cj f/3 and S(^) 0, sec Fig. |. The coefficient 
of effective uniaxial anisotropy A{£,) slowly depends on ^, 
namely yi(0) = (Vf + 4£2-l)/(2£)-2/3 and ^1(1 ) = 1/3. 
When the particle aspect ratio £ < f then A{£^) < 0, see 
Fig. ^ and we have an effective easy-plane anisotropy. 
When £ > f, then A{£,) > and we have an effective 



easy-axis anisotropy. More details are given in Sec. [II A. 

In addition to the effective uniaxial anisotropy given 
by A{(,), we have the essential 23 (^) term which gives an 
effective in-plane anisotropy. For the disk-shaped parti- 
cle this anisotropy coefficient is always real, arg23 = 0. 
The value of 23 is almost at origin but its contribution 
becomes important at the boundary, see Fig. ^. We ob- 
tain the following asymptotics, valid form small e and 
1/2<?<1 



arctan ( - 



2£(3e - 2) 



In 



16 



1-e 
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£2 + (I-0' 



£2 + (I-e)2 

(25) 



(see Appendix Thus the 23 (^) term causes bound- 
ary effects and is responsible for the configurational 
anisotropy. In the limit £ ^ (more precisely, when 
ao h <^ R) the 23 (^) term is concentrated near the 



boundary, corresponding to the surface anisotropy. 

The energy of the nanodisk can be derived from 
Eq. (|l|): 



= Wex + Wd, 

R 

■Itt 1 



(26) 



R 'Itt 

rdr J dx 









In the next subsections we analyze the homogeneous state 
and the vortex state. 



A. Homogeneous state 

Let us consider a homogeneous magnetization along 
the X direction of the disk-shaped nanodot, so that 9 = 
7r/2, (/) = 0. The exchange energy vanishes. The second 
term in the dipolar energy (|2^) also vanishes because of 
averaging on x- The total energy is then 



1 

= 2^jAmd^ = W^s{e) 





27r 



W^Ss(e) = ^{-l + yT+^[£'K(m) 



l-£^ E 



(27) 



where to = (1 + £^)^^, K(?ti) and E(to) are the com- 
plete elliptic-integrals of the first and the second kind, 
respectivelycj. The constant term — 27r/3 is the isotropic 
contribution. The second term W^g is the well-known 
magnctostatic energy of the homogeneously magnetized 
disk, first calculated by JosephE£l. 
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Disk diameter cross section, (r/l) 

FIG. 3: (Color online) Comparison of the vortex profiles 

for the micromagnctic simulation and the effective 
anisotropy approximation for a Py nanodisk (2i? = 212 
nm, h = 16 nm). The red curve corresponds to the 
spin-lattice simulations for the effective anisotropy 

model with ?{ = J{ex + ^d"- The blue curve 
corresponds to the micromagnctic simulations. The 
black dashed curve to the gaussian ansatz 
cos 6* = exp(— r^/r.y). 



If the disk is now homogeneously magnetized along the 
z-axis, then = 0. From (^6|) one sees that the corre- 
sponding total energy = —2W^. The transition be- 
tween these two homogeneous ground states occurs when 
= W^. This happens only for Wx = 0, i.e. for 
^Ms(^c) = 27r/3. This gives a critical value Sc ~ 0.906 
which agrees with the result by AharoniEH. 



B. Vortex state 

Let us consider a nonhomogeneous state of the disk- 
shaped particle. In this state the system has a larger ex- 
change energy compared to the homogeneous state. This 
should be compensated by the dipolar term. According 
to the dipolar interaction always favors a spin dis- 
tribution of the form 



X± 



(28) 



where we take into account that the in-plane anisotropy 
constant 23 takes real values only. Such a configuration is 
called a vortex. In highly anisotropic maeaets there can 
exist pure planar vortices with 6 = 7r/2B3 However we 
consider here out-of-plane vortices, realized in "soft" ma- 
terials typical of nanodisks. The out-of-plane component 
of the magnetization has a radial symmetric shape, and 
it almost does not depend on z for thin disks, 9 = 9{r). 
Now we can calculate the vortex energy. The vortex solu- 
tion ( p8| ) is characterized by cos 2(0— x) = — 1, providing 
the minimum of the in-plane component of the dipolar 



energy: 



VF'^-27r/^dC 







The exchange energy term 

2 



Wl 



in 
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rdr 





Finally, the vortex energy is 

rvortcx jirx , jr/vortex 



sm 



2 a-\ 



.(29) 



(30) 
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1/1/ EP 
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EP 
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sin^ 9 cos^ 9 



F(e) 27r ^ Cd^j [3^1(0 + 2] cos^ 9(R^) 
+ 23(0 sin2 0(i?e)}. 



(31) 



Avortcx 
EP 

an easy-plane magnetJ 



Here W^^^'^ coincides-iwith the energy of the vortex in 



vortex 
EP 



2^( - 

R 



In 



A = 5.27 (32) 



and Fie) is the configurational anisotropy term. The vor- 
tex state is energetically preferable to the homogeneous 
state when the configurational anisotropy term exceeds 
the energy of the easy-plane vortex F(e) > W^^^^. This 
relation allows to calculate the critical radius i?c by solv- 
ing the equation 



27r 



- 1 In 

R 



tKR^ 



= Fie). 



(33) 



To calculate the integral in F{e) we use the trial function 
for the vortex structure 



= cos 6* = cxp(— r^/r^). 
The core width depends on the disk thickncss0 
r^{h) « £V2^l + ch/£, c w 0.39. 



(34) 



(35) 



The relation ( |33| ) providing the border between the easy- 
plane and the out-of-planc vortex states can be analyzed 

in the limit e 0. Then F{e) - (27r£/3) ln(^7r/(2e)) , 

hence R^"^ w i^/S/l^^his is in qualitative agreement 
with previous resultsElrEil. 

Let us estimate now the contribution of the dispersive 
part of the dipolar energy. Taking into account that for 
the curling state (p|) the second term in Eq. ( 23 ) vanishes: 



V • m = and that the Fourier-component 



of the out-of-plane component (^ 
TTrl e~'?^''t', from Eq.(^) we get 



has the form m~ 




for ry > h, 
for ry <C h. 



(36) 
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(Col or o nline) Numerical results for the vortex state P y pr ism (sides 212 x 212 nm, thickness /i = 16 nm 



a] and |(b)| represent the spin-field distribution, and Fig. 



the configurational anisotropy lines. These lines 
determine the in-plane anisotropy axis orientation, calculated from Eq. ([2l|); lines lengths correspond to the 

anisotropy amplitude in a particular point. 



Comparing Eq. ( p3q ) with Eqs. (31) and (p2[), and taking 
into account (^5[), one can conclude that in the limit e 
the dispersive part of the dipolar interaction docs not 
change significantly the vortex stability criterion. More 
precisely, our effective anisotropy approximation works 
correctly not only for e ^ but also for h, i? ^ 
Our numerical results show that it gives the vortex state 
as an energy minimum for disk diameters 2R > 30^. 



IV. NUMERICAL SIMULATIONS 

To check our effective anisotropy approximation, we 
performed numerical simulations. We used the pub- 
licly available three-dimensional OOMMF micromag- 
netic simulator codc£3. In all micromagnetic simula- 
tions we used the following material parameters for Py: 
A = 1.3 X 10-s erg/cm (using SI units A^^ = 1.3 x 10"" 
J/m), Ms = 8.6 x 10^ G (Mf = 8.6 x 10^ A/m), the 
damping coefficient rj = 0.006, and the anisotropy has 
been neglected. This corresponds to an exch ange length 
£ = y^A/inM^ « 5.3nm {(^^ = y^A/fXoM^). The mesh 
cells were cubic (2 nm). 

We also test our effective anisotropy approach by the 
original discrete spin-lattice simulator. The spin dynam- 
ics is described by the discrete version of the Landau- 
Lifshitz equations with Gilbert damping 



d5„ 

dt 



d'K 



dSr, 



1 
S 



dt 



(37) 



which we consider on a 2D square lattice of size (2i?)^. 
We have assumed a plane-parallel spin distribution ho- 
mogeneous along the z-direction. Each lattice is bounded 



by a circle of radius R on which the spins are free cor- 
responding to a Neuman boundary condition in the con- 
tinuum limit. We integrate the discrete Landau-Lifshitz 
equations ( p7|) with the Hamiltonian 5f = K^x+^d given 
by (H) and ([7a|), using a 4th-order Runge-Kutta scheme 
with time step O.Ol/TV^. These spin-lattice simulations 
were done to validate our analytical calculations for the 
effective anisotropy model. Throughout this work we 
compared the results of the spin-lattice simulations with 
IK = IKex + with the results of micromagnetic sim- 
ulations. We never found any noticeable difference. We 
present the results for a disk-shaped and a prism-shaped 
nanoparticle because these two geometries are the most 
common ones in experiments. 



A. Disk— shape nanoparticle 

Our effective anisotropy approximation provides the 
exact solution for all homogeneous states for a nanodisk. 
Therefore we do not need to justify it for the homoge- 
neous states. We consider here the vortex state. As 
we analyzed before, the model can provide the prefer- 
able vortex state for disk diameters 2R > 30£, which is 
in an agreement of the model usage critcrium (^. We 
compare the magnetization distribution in the vortex for 
our effective anisotropy model and for the the micromag- 
netic simulations. Since the in-plane vortex structure is 
characterized by the same distribution (jj = x ± 7r/2 for 
both methods, we are interested in the out-of-plane vor- 
tex profiles. We performed such a comparison for a disk 
of size 2R/£ — 40 and h/£ — 3, which satisfy the cri- 
terium (||). The results are presented on Fig. One can 
see that the vortex shape from the effective anisotropy 
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;r 3;r 2ff iL ;r 3^ 2;r 

2 2 2 2 

X X 



(a) Distance from the vortex center (b) Distance from the vortex center 
r = 10£. r = 201. 

FIG. 5: (Color online) The in-plane spin angle as a function of the polar angle for the vortex state in a prism of 
Py of sides 212 x 212 nm and thickness ft. = 16 nm. The red dashed curves correspond to the effective anisotropy 
approximation, the blue solid curves — to the micromagnetic simulations data. 



model agrees with the one obtained from the micromag- 
netic simulations within 0.11 in absolute error. 



V. DISCUSSION 



B. Prism— shape nanoparticle 

Now we check the validity of the effective anisotropy 
approximation for the prism-shaped nanoparticle. We 
chose this shape because there are numerous experiment 
with a square geometry, see for a review Ref. |l} We per- 
formed the two types of simulations for a square shaped 
nanoparticle, see Figs. ^ and The two equilibrium 
magnetization distributions, obtained for the micromag- 
netic model and the spin-lattice simulation agree with a 
very high precision. 

As discussed above the large scale distribution of the 
magnetization is described by Eq. (^). Calculating 
numerically the coefficient 23 (see Appendix ^ for de- 
tails), we found the distribution of the configurational 
anisotropy lines for the square geometry. This is shown 
in Fig. ^ The comparison of Figs. |^ shows that the ef- 
fective anisotropy lines corresponds to the magnetization 
direction in the main part of the system. Note that the 
effective anisotropy approach fails near the corners: the 
sharp field distribution near the prism vertices (Fig. ^) 
is not energetically preferable when the exchange contri- 
bution is taken into account. 

We can also check the validity of the effective 
anisotropy approach for the complicated "vortex" struc- 
ture in the square geometry, by comparing the distribu- 
tion of the in-plane spin angle (j) to the one given by the 
micromagnetic simulations. This is done in Fig The 
figure shows that the two different approaches agree very 
well. The (/)(x)-depcndencics coincide within 0.11 in ab- 
solute error for r = 10^ and within 0.04 for r ~ 201. 



To summarize, assuming that magnetization is inde- 
pendant of the thickness variable z, we have reduced the 
magnetic energy of a thin nanodot to a local 2D inho- 
mogeneous anisotropy. The first term A determines the 
uniaxial anisotropy along the z-axis. The second term 23 
gives the anisotropy in the XY-j>\aiie. 

For thin nanoparticles e < 1 the term A w const < 0, 
gives an effective easy-plane anisotropy. This general- 
izes the rigorous results obtained for infinitesimally thin 
filmsEj. The function 23 (x, y) is localized near the edge of 
the particle so that spins will be tangent to the bountiapj;. 
This confirms the notion of a surface edge anisotropyll3ll3. 
When the nanoparticle is thick £ ^ 1, the anisotropy con- 
stant yi > 0, is again almost constant and the spins will 
tend to follow the z axis (easy-axis anisotropy). The in- 
plane anisotropy 23 depends on the thickness, see Fig. |2^. 
The special distribution of 23(a;,y) is responsible for the 
volume contribution of the dipolar energy. 

The above effective anisotropy approach: (i) shows the 
nature of the effective easy-plane anisotropy and the sur- 
face anisotropy, (ii) generalizes the surface anisotropy for 
the finite thickness, and (iii) gives a unified approach to 
study dipolar effects in pure 2D systems and 3D magnets 
of finite thickness. 

It is instructive to make a link between our approach 
and the rigorous results which were obtained in Refs. [T^ , 
|l|,|l3jl8[|l|,|^,|l[ Our equations (|3^), (|3|) show that for 
the vortex ground state to exist, it is crucial to have both 
types of anisotropy: out-of-plane anisotropy and in-plane 
one. It is shown by Kohn and SlastikovEj that the energy 
of a thin magnetic film with an accuracy up to can be 



9 



presented as the sum 



E = £^cxch 



^bdry -^trans 

-21 In 



27r 



2/2 (^^) 

(m • n) + £ / (w^) . 



Considering the Umit e ^ and ^^/(el lne|) = const we 



see from Eq. (^8|) that formaUy the last term is dom- 
inating and its contribution has to be accounted as a 
constraint = 0, see Ref. |l9[ This constraint prevents 
the existence of the vortex ground state of the nanodot 
because the energy of the vortex in the continuum hmit 
is infinite due to divergence at ?■ ^ 0. However this di- 
vergence is removed by the out-of-planc component of 
vortex which is described by a locahzed function with 
radius of locahzation ^ i [see Eq. (35)]. This means 
that the last term i?trans in (|8|) scales like the exchange 
term ii^cxch- In this limit all three terms of (^8|) are of 
the same order and provide the existence of the vortex 
ground state. 

This reduction of the nonlocal dipolar interaction to 
a local form is a first step towards an analytical study 
of nanomagnctism. Wc developed a method of effective 



anisotropy and illustrated it on a few examples. We plan 
to apply this method to the dynamics of vortices in nano- 
magnets. 
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APPENDIX A: DISCRETE DIPOLAR ENERGY CALCULATIONS 



Let us consider the dipolar interaction term Vi^. Using the notations 



^nm ynm -^nm 



ao ao 
one can rewrite this energy as follows: 



ao 



o2 + 22 



(Al) 



Z QZ ,2 

Tim 



n.m \ 



r3 

' nm 



' nm 



{S^Xnm + S^Unm) {S^Xnm + S^Unm) f ^ ^ \ St^li^ z{pv fj.) + {S^S^ + S^S^^ Kl{pvfj.) (A2) 



{S^Xufj, + S^iji,^) {S^Xufj, -\- S^Hi^n) K2{pvfj.) 



^vfi^ Unm, 



Here we used the obvious relations Xnm 
not depend on the z-coordinatc: Sn = S^,, S. 
kernels Ki, K2 and contain information about the original 3D structure of our system. 



VuiJ. and the basic assumption that the magnetization does 
S^. This allows us to reduce the summation to the 2D lattice. The 



1 



1 



2 X! , ,2 ^3/2' 



(s^ + ^nm) 



K2{s) 



1 



1 



2zl 



^ (,2+;2 )5/2' ^^(^) = 2 ^ (,2 + ,"2 "7/2 ■ (A3) 
nz,7nz \ ' nmJ 71^,771^ ' '^nmJ 



Taking into account that 
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one can present the dipolar energy in more symmetrical way: 



(A4) 



The total Hamiltonian is the sum of two terms (I) and (0). 

Here we show that the main effect of the nonlocal dipolar interaction is an effective nonhomogeneous anisotropy. 
Using equality 



where Cnm — Cmm one can split the dipolar Hamiltonian (A4) into a local contribution and a nonlocal correction 



^K2{p.^)x.^y.^ {sz - s^^) {sy - sy^] 



(A5) 



APPENDIX B: CONTINUUM LIMIT OF THE LOCAL DIPOLAR ENERGY 

Here we present the continuum limit of the discrete dipolar Hamiltonian (|^) corresponding to the dipolar energy 



(Bl) 



where m^, ~ '^if^v- Hence the continuous magnetization vector m according to Eq. (^2|) takes the form m(r) 
X^i^ »T^i^<5 (r — Ti,). Here A^, Bi, and are determined as follows 



.2 



^ g ^ ^ 2;'^ y"^ 



C^^ Y K2{p„^)2x^ij,y^fj, = ^ E E 



2xjj^yi/f^ 



2 V ^ f«2 , ^2 v5/2- 



(B2a) 

(B2b) 
(B2c) 



The continuum version of the effective anisotropy constants ( |B2| ) can be found using a relation 

E E^(l^"-l)^-/ Az'F{\z-z'\) + - Az[F{\z\) + F{\h-z\)]+-[F{Q) + F{\h\)] 



n, — m~ —0 



4 / dzi^(|z|) [/i - z + ao] + i [F(0) + F[\h\)] , h = N,ao > 0. 
o-o Jo ^ 



(B3) 
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Let us start with the calculation of the coefficient A^, from Eq. ( B2a ) 
A{x,y) 



at 

2TTh' 



A,, 



1 (Ai +A2+ As) , Ai = lim 



d^x' I dz 





2 /•2ir 



da ao _ al 



d^x 



\r—r' I >r * 



(p2 + z2)5/2 

2 '•^^ PMa 



(B4) 



(p2 + /j2)5/2 8^7o (P2+/l2)3/2 4(^2 + /,2) 3/2 



Here p = •^/(a; — a;')2 -|- (y _ y')2 g^j^^^ -^g used a local reference frame (^6|) and the Heaviside function 0+(x) takes the 
unit values for any positive x and zero values for x < 0. The Heaviside function is added here to fulfil the condition 
yii = in a 2D case, when for h = 0. There is a singularity in Ai, due to the nonintegrability of the kernel at 
^nm = 0. To regularize it we use a method similar to the one in Ref. Specifically, we present Ai in the form 
Al ^ Al — Aq. The coefficient Ai is a regular one: 



Al = 



, .\ i2z^-p^){h-z + ao) n 1 

d X I dz — = ^ aQ<a+(h) 





271- J Jo " (p2+z2)5/2 

The singularity is inside the ^Iq term: 



271 



da 



VW+h^ - P 



aoh 



Ao = 



Q+{h) 



lim 



27r r*^0 

|r— r |<r* 
z=0, z'>0 



d-'x'd 



,^ ,(2z'' -p'){h^z' + aa) e+ih) 



n2 _i_ ^/2^j5/2 



Ip^ + z' 



2tt 



(h + ao)/i - h 



Ii = lim 

r*^0 



2j-'2 - rt2 
j2 /j„/ ^2 p 



d^x'dz' 



\r — r' I <r* 
z=0. z'>0 



(p2 + z'2) 



lim / d^x'dz' ^ 



5/2 

|r — r'|< 



lim 



dz'^ \r - r'l 3 »-*^o 



d^x'A-. 



r — r' 



z=0, z'>0 



\r — r' I <r* 
z=0, z'>0 



linr / d^x'5{r-r') = - — , 

r — I <r* 
z=0, z'>0 



/t = lim 



d2^'dz'f^i?f^ ^ = — lim / d3x'z'5(r - r') = 0. 

5/2 3 r*^0 ' ^ ' 



r— r |<r 
z=0, z'>0 



(p2 + z'2) = 



\r — r' I <r* 



(B5) 



Finally, Aq = — [/i + ao0+(/i)] /3 and the coefficient of effective anisotropy A{x,y) takes a form (15a). 



The coefficients B^, and C^, can be calculated in the same way, starting from Eq. (B2b) 



4p2ix 



age 



2'Kh 



3ao 
47r/i 



\5/2 /i 







dae 



(B6) 



-2{h + ao) In 



-2lQ 



2 /•27r ^-2i, 



3a^ 



'da 



^ / dae- 



2 r2n 



o 2 rt2p-2'" 
StT y (p2 + ;j2)5/2 StT 7o 



VP2 + fe2 _ 
P 

3P2 + 2/i2 



2jq 



(P2 + /i2)3/2- 



Finally, the coefficient of effective anisotropy ^(x,?/) takes a form ( |l5b| ). As a result the dipolar energy (Bl) can be 
expressed as (^. 

Note that for the circular system one can obtain exact expressions for the coefficients A and 23. Let us first find 
the coefficient A. Assuming that ao (or equivalently ag — > 0), one can rewrite the coefficient A, see Eq. (B4), as 
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follows: 



Ia{x) 



1 1 

3 Aire 
2 



da 



lA{2e)-lA{0) , Ia{x)- , , ^ , , 

V^^ + e' +.i;2-2ee COS a 

=_{ [x^ + + 1)'] E(m) +[l-x^- e] K(m) + + F^{x)\ - 27rx, 

C + 1) 

4^ 



(B7) 



x^ + {i + 0'- 



2^ 



where n(z/±|/i) is the complete elliptic integral of the third kind.E3. 

To calculate the in-plane anisotropy coefficient 23, see Eq. (|B6|), it is convenient to use the following relations 



Re [Be-^^^] = -^B^ = / - z)I.{x), h{x) = 3 / d 



, {x-x'f-{v-y'f 



(B8) 



92 92 



1 



^dydy' dxdx' J ^ + 
For a circular system dZ' = i?dx' (— e^jSin^' + ej,cosx'), hence 

9 . , 5 



(/92 + ;22)5/2 

/ [V X F] • do- = / F F = X V , 



^(x - x'Y + {y- y'Y + 



hix)=Rjdx' 



TP sm X - ^ cos X 
ay ox 



1 d 

= ri?cos(2x) — 

^r2 + i?2 _ 2ri?cos(x - x') + ^2 



27r 



COS a da 



?'7 ^7.2 + i?2 - 2ri? cos a + z2 



Taking into account that ImS = for the circular system, one can calculate finally the effective in-plane anisotropy 
coefficient S as follows: 



1 

27r£ 



lB{2e) - Ib{0)\ , Ib{x) = ciK(/^) + C2E(^) + cgH 
2 



4e 



(1 + 0^ 



2-2x2-e2_ ^.2^^2)- U2 + ^2 _ 2) ^^2 + (1 + 
Cl = , ^ , C2 = ^ ^— T5 , C3 



(B9) 



ie^/x^T{T+W ' ' C2(l+^)^^2 + (^+l)2- 

The dipolar energy Wd [see Eq. (p6|)] for the disk-shaped system can be presented in the form Wa = + Wd, where 



Wd = -^ I d^x 



A{r)+'B{r) cos2((/.-x) 



sin^ei 



(BIO) 



and Wjl = —2R ^ J d'^xA{r) being the isotropic part, the effective easy-plane anisotropy parameter A — 3A. 



APPENDIX C: CONFIGURATIONAL 
ANISOTROPY FOR A HALF-PLANE AND A 
SQUARE PRISM 

We start here with the problem for a half-plane. Con- 
sider the large scale behavior of the dipolar energy, 
given by the in-plane effective anisotropy 'B(x,y), see 
Eq. (|2^). Straightforward calculations lead to the effec- 
tive anisotropy constant for the upper half-plane 



1 r°° 

■B{x,y) = S(2/o) = ^J dxyoiyl " 



pi 



2TTh 



In 



yl + h^ 



1 h 
— arctan — , 
71" Vo 



(CI) 



where we choose the origin of the local reference frame at 
the boundary of the domain, at (a:, y) = (0, 0), yo denotes 



the distance from the boundary, and P = \/ x^ + 2/o ■ ^^^^ 
can see that S does not depend on x, it takes only pos- 
itive real values, hence argS = for any distances y^ 
from the boundary. This means that the in-plane spin 
angle is al way s para llel to the half-plane edge. Using 
Eqs. (20b), (20c) and (CI) we found that the main con- 



tribution to (CI) is provided by the boundary domain 

/yoTi. Since this domain col- 
't 0, we conclude that for any 



X € [— i?o; ^o] with i?o ~ 
lapses to a point when ?/o 
geometry the in-plane spin distribution is parallel to the 
boundary near the edge. If the curvature radius of the 
sample boundary is larger than i?Oj then spins are par- 
allel to the boundary over a distance smaller than i?g//i. 
One should remember, that this conclusion is adequate 
for regions, where exchange interaction has no principal 
influence. 
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FIG. 6: (Color online) Arrangement of coordinates in 
the local reference frame for the prism shaped particle. 



Let us consider now the configurational anisotropy for 
the square prism, which has the diagonal 2R, see Fig. ^. 
It is convenient to use the local reference frame in the 



same way as in Sec. II A. The relative polar coordinates 
arc defined as follows: 



Rn ^R^l+e ~2^cos{mT/2-x), 



(fn = arccos — 

RnRn+l 



P, 



+ Rn+1 

2RnRn+l 

sin Lp„ 



where ^ = + /R- Now we are able to compute 
magnetization distribution on a large scale, which follows 
from the minimization condition (^ij). Straightforward 
calculations give 



X + 2 - 2 ' 



S — 



1 

2ti 



Jipo-Va 



3 

j=i •'i'l-i 
j 



37r ■ f Ro sin tpo 

Y — arcsm r= — 

4 V RV2 



(C2) where 3^iPi,h) is defined by (p03) 



(C3) 



(C4) 

(C5) 
(C6) 



i?V2 cos (q + x - (2jt+ l)'/i/4)' 
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